Topological Signature of First Order Phase Transitions 



(N 
O 

o 



43 

cd 

CD 



o 

cd 



> 



o 



- 1—1 

X 



L. Angelani 1 ' 4 , L. Casetti 2 , M. Pettini 3 , G. Ruocco 1 , F. Zamponi 1 . 

1 Dipartimento di Fisica and INFM, Universita di Roma La Sapienza, P. A. Moro 2, 00185 Roma, Italy 
2 INFM - UdR Firenze, via G. Sansone 1, 50019 Sesto Fiorentino, Italy 
3 Istituto Nazionale di Astrofisica, Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 1-50125 Firenze, Italy 
4 SMC-INFM, Universita di Roma La Sapienza, P. A. Moro 2, 00185 Roma, Italy 

(Dated: February 6, 2008) 



We show that the presence and the location of first order phase transitions in a thermodynamic 
system can be deduced by the study of the topology of the potential energy function, V(q), without 
introducing any thermodynamic measure. In particular, we present the thermodynamics of an 
analytically solvable mean-field model with a fe-body interaction which -depending on the value of 
£f>l . k- displays no transition (fc=l), second order (fc=2) or first order (fc>2) phase transition. This 

rich behavior is quantitatively retrieved by the investigation of a topological invariant, the Euler 
characteristic x(v), of some submanifolds of the configuration space. Finally, we conjecture a direct 
link between x( v ) an d the thermodynamic entropy. 



PACS numbers: 05.70.Fh; 02.40.-k; 75.10.Hk 



Despite their major physical relevance, first-order phase transitions are still lacking a satisfactory theoretical un- 
derstanding of their origin. More specifically, nothing comparable to the renormalization group analysis for critical 
phenomena exists in the case of first order transitions, and, on a more mathematically rigorous background, neither 
in the Yang-Lee theory for the grand-canonical ensemble Q nor in the Ruelle, Sinai, Pirogov theory in the canonical 
ensemble [0 an a priori mathematical distinction can be made among the interaction potentials leading to first or 
second order phase transitions respectively. 

The present work aims at contributing a step forward on this topic. By simply knowing the microscopic interaction 
potential, thus even prior to the choice of a statistical measure, we seek a characterization of the occurrence of a phase 
transition and of its order. As we shall see in what follows, this can be actually achieved within the recently proposed 
topological approach to the study of Phase Transitions (PT) ||, [|, || . Within this approach, that at present has dealt 
only with second order PTs, the existence of a PT is signaled by a singular point in the energy dependence of some 
topological invariant. In particular, in the specific cases of the mean-field XY [Q and of the ip 4 models it has been 
shown that it exists a discontinuity in the first energy derivative of the Euler characteristic, x( u ): of some submanifolds 
, of the configuration space, located at that energy, v c =v(E c ), where the system undergoes a second order PT. The 
qq ■ quoted works demonstrated, therefore, that important thermodynamic features -the second order PTs- are strictly 
related to the topology of the potential energy hyper-surfaces. In particular the knowledge of the energy distribution 
■ of the critical points (saddle-points) of the potential energy function V(q) allows one to predict the existence and 
CO ' location of second order PTs. 

The -still lacking- extension of the topological approach to first order PTs is obviously of great interest. It would 
allow to get new insight into the challenging problem of the origin of first order PTs, providing a new method 
of tackling them, and it would validate the topological approach as a possible theoretical method of unifying the 
treatment of different kinds of PTs. The latter point is particularly interesting in view of encompassing also more 
"exotic" transitions, among which it is worth to mention the glass transitions. In this respect, we notice that indication 
i-^j ■ on the existence of a relation between the glass transition and the underlying topology of the potential energy surface 
already exist. Indeed, it has been shown || that, in simulated realistic glasses, the order of the saddles visited at a 
O ■ given temperature T vanishes when T approaches the (dynamical) transition temperature. 

CD ■ l n this work, we show the existence of a relationship between the topology and the thermodynamics of a dynamical 
system. In particular, we introduce and discuss here the " fc-trigonometric model", an analytically solvable mean-field 
model with a fc-body interaction. According to the value of the parameter k, the system has no PT (fc=l), undergoes 
a second order PT (fc=2) or a first order one (fe>2). Remarkably, for this model an exact analytical computation 
is possible of both the thermodynamics and of a fundamental topological invariant, the Euler characteristic \(v) of 
the submanifolds M v of the configurational space defined by V(q)<Nv. In the topological properties, the existence 
of the PT is signaled, for fc>l, in the usual way [Q, ||: a discontinuity of the first derivative of xi v ) is observed at 
v c =v(E c ), i. e. at the potential energy level where the PT actually takes place. Remarkably, this happens also in 
the case of the first order PT. But for k — 2 the second derivative d 2 x(v)/dv 2 is negative around the singular point, 
while for k>2 it is positive. Therefore: i) A singularity in the first derivative of x(v) is a signature of the PT, both 
in the first or second order case, ii) The sign of the second derivative of x( v ) around the singularity allows to predict 
the order of the transition. The latter conclusion leads to a further, important implication: it is well known [7| that 
regions with d 2 S(E) / dE 2 > are not thermodynamically stable (they would correspond to negative specific heat) 



2 



and their presence corresponds to the existence of a first order PT; therefore, we conjecture a relation between the 
thermodynamic entropy and the topological Euler characteristic. 

We study the properties of a specific model, the mean field k- Trigonometric Model (fcTM), defined by the Hamil- 
tonian 



H k (p, q) = K{p) + V k (q) =J2^r + ^*(«) ' W 



2m 

I — 1 

where the potential energy is given by 

v k(q) = j^T [i-cos^(q n + ... + «<»)] ; (2) 

A and L are the energy and length scales. If k=l the fcTM reduces to the Trigonometric Model (TM) introduced 
by Madan and Keyes | as a model for the Potential Energy Surface (PES) of simple glass-forming liquids that 
reproduces many properties of the "true" PES of these systems Q. For a given k, the model is a mean-field model 
with fc-body interaction. It is convenient to introduce the angular variables <pi=(2ir / L)qi. The model has then a 
symmetry group C kv obtained by the transformations ipi — > ipi + 2irl/k and ipi —> —ipi. If we think at ipi as the 
angle between a unitary vector in a plane and the horizontal axis of this plane, these transformations are rotations 
in this plane of an angle 2irl/k and the reflection with respect to the horizontal axis respectively. Using the relation 
cos((fi 1 + ... + (fii k )=Re(e tVi i ■ ... ■ e 4Vi fc), the potential can also be written as: 

V k (<p) = NA[1 - Re(c(p) + is(ip)) k ] , (3) 

where c{tp)=jj cos <pi and s((p)=jr sin^;. As usual in simple mean field models it is easy to calculate the 
microcanonical partition function, given by 

r d N v d N a 

^N,k{E) = J ±^±Ss(H k -E). (4) 
Using the integral representation of the delta function, we get 

« = /^T e "' M - (5) 

Now, as we are looking for a saddle-point evaluation of the integral over /?, we can rotate the integration path on the 
imaginary axis in the complex-/? plane. This is justified because, as we will check at the end, the saddle-point is located 
on this axis We can now perform the integration over the momenta and use the fact that V k (ip)—V k (c((p),s(if)), 
see Eq. (||), to obtain 

n Ntk (E) = C N P N J dp di d v e P(E-v k (S, v )) . 

■Jd N <pS[N(Z-c(<p))} 5[N(ji-8(<p))] , (6) 

where p=N / L and the constant Cjv gives only a constant contribution to the entropy per particle, i. e. it is at most 
of order e N . The last integral can be evaluated using again the integral representation of the delta function, and 
rotating then the integration path as previously discussed; it turns out to be: 

dfJ, dv c _ N{lli+vri) I dN ^^cosipi+umji^i) = 



(2nf 



^e- N ^+^(2nI (A)f , (7) 



(2tt 



having defined K=y/ ji 2 + v 2 and the Bessel function 



We can then write the partition function as 

n N , k (e)=C N p N I dme Nf ^ m ' e \ (9) 
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where m=(/3, £, T], fx, v), e—E/N and 

f k (m, e) =j3e - /?A[1 - Re(£ + «7) fc ] 

1 (10) 
- - log/3 - /if - vt) + log Jo (A) . 

Then, using the saddle-point theorem, the entropy per particle, s—S/N, is given by (fcs=l): 

Sfe(e) = lim — logfi;v,fc(e) = max/t(m, e) . (11) 

N^oo iV m 

To find the maximum of /fe(m, e) one can calculate analytically some derivatives of / to obtain a one-dimensional 
problem that can be easily solved numerically. In Fig. 1 we report the caloric curve, i. e. the temperature vs energy 
relation T(e)=[ds/de]~ 1 for three values of k, k=l, 2 and 3. As expected, the temperature is an analytic function of 
e for k=X, because in this case there is no interaction between the degrees of freedom. For k—2 the system undergoes 
a second order phase transition at a certain energy value e c , that changes to first order for k>2. At the transition 
point the Ckv symmetry is broken, the order parameter being the "magnetization" m=(\c + is\); for e<e c there are 
k pure states related by the symmetry group. It is interesting to note that -for the present model- this happens also 
in the case of the first order transition. 

Having an exact solution of the thermodynamics of the system, we can look for some topological signature of the 
phase transition. In ji], || it was stated that the Euler characteristic x( v ) of the submanifolds M v ={q \ V(q) < Nv} 
of the configurational space shows a singularity in correspondence of the potential energy value v c =v(e c ) at which 
the transition takes place. The result of || was obtained in the case of a second order PT; our aim is to extend 
this result to the case of a first order PT. Remarkably, the Euler characteristic of M v can be calculated analytically 
in our model. The general definition is flj0| : 

N 

X (v) = X (M V ) = ^(-1)>„(M„) , (12) 

n=0 

where the Morse indexes jj, n (M v ) correspond to the number of critical points of order n of the function V(q) that 
belongs to the manifold M v . The critical points (saddles) (p are defined by the condition dVk(<p)=0, and their order n 
is defined as the number of negative eigenvalues of the Hessian matrix TC^(ip)=(d 2 Vk/ difid(pj)\0. To determine the 
location of the saddles we have to solve the system 

— - = -A k Hje[i(c + is) k - 1 e i ^} = 
dip, 

A k C fc_1 smp - 1)V> + tfj] = , Vj , (13) 

where we defined c+is—^e 1 ^. From Eq. (||) we have Vjt(<^)=iVA[l — £ fe cos(kip)]; then the saddles with £(y?)=0 have 
ener gy y =V(0)/N=A. We can neglect them because, as we will see at the end, v—A is a singular point of x( v )- Then 
Eq. (|13j ) becomes 

sin[(fc - l)xl> + if,] = , Vj , (14) 

and its solutions are 

Cpj = [UjlT ~(k- l)V']mod 2tt , (15) 

where njS{0, 1}. The saddle point tp— is then characterized by the set n={rij}. To determine the unknown constant 
ip we have to substitute Eq. (JL5|) in the self-consistency equation 



= c + is = N- 1 e ' Vl = N- l e-**^-^ • ( i6 ) 

3 3 

If we introduce the quantity x(<p) defined by 

x = N~ 1 J2 n 3 1 - 2x = A r_1 ' 

3 3 
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FIG. 1: Microcanonical temperature as a function of energy for three different values of k; for k=l there is no phase transition, 
while for k=2 there is a second order transition and for k>2 a first order one. In the inset, the potential energy v vs total 
energy e curve; the critical point is, for Vfc>2, v c =A. 



we have from Eq. (\i 



C = 1 1 — 2x1 



4>i 



2ln/k 

(2/ + l)vr/fc 



for x < 1/2 , 
for x> 1/2 , 



(18) 
(19) 



where I GZ. Then the choice of the set {fij} is not sufficient to specify the set {<fij} because the constant ip can 
assume some different values. This fact is connected with the symmetry structure of the potential energy surface (the 
different values of ipi generate the multiplets of saddles). We have then obtained that all the saddles of energy u^A 
have the form 



fif' = i n J n — ( k — l)V'/]mod 2tt ■ 

A k Re[N- 1 (k - l)(c+ i s ) k - 2 e i ^ + fi 



The Hessian matrix is given by 



In the thermodynamic limit it becomes diagonal 



U% = Sij A k C^ 1 cos (ip(k - 1) + ipi) 



(20) 



(21) 



(22) 



One can not a priori neglect the contribution of the off-diagonal terms to the eigenvalues of 7i, but we have numerically 
checked that their contribution change at most the sign of only one eigenvalue over N. Neglecting the off-diagonal 
contributions, the eigenvalues of the Hessian calculated in the saddle-point tp are obtained substituting Eq. (|2^) in 
Eq. (H): 



A, 



(-l)^A k C 



k-l 



(23) 



so the saddle order is simply the number of %=1 in the set n; we can identify the quantity x{(p) given by Eq. ( |17| ) 
with the fractional order n/N of (p. Then, from Eq. (||), ( |l8| ) and (jlj]) we get a relation between the fractional order 
x((p) and the potential energy v(<p)=V(ip)/N at each saddle point (p: 



x(v) 



1 — sgn 1 1 



jn(l 



v 
A 



i/fc' 



(24) 



v/A 



FIG. 2: a(v) (see text) as a function of the potential energy. The phase transition is signaled as a singularity of the first 
derivative at v c =A; the sign of the second derivative around the singular point allows to predict the order of the transition. 
The region v>A, in which a'(v)<0, in not reached by the system (see the inset of Fig. 1). 

Moreover, the number of saddles of given order n is simply the number of way in which one can choose n times 1 
among the {rij}, see Eq. (|o|), multiplied for a constant Ak that takes into account the degeneracy introduced by Eq. 
(|l9|). Therefore: i) the fractional order x—n/N of the saddles is a well defined monotonic function of their potential 
energy v, given by Eq. (pi]), and ii) the number of saddles of a given order n is Ak( N )- Then the Morse indexes 
fi„(M v ) of the manifold M v are given by Ak^) if n/N<x(v) and otherwise, and the Euler characteristic is 



It has to be stressed that o~(v) is a purely topological quantity, being related only to the properties of the potential 
energy surface defined by 14 ((/?), and, in particular, to the energy distribution of its saddle points. From Fig. 2 we 
see that there is a signature of the PT that is evident in the analytic properties of <r(v). First, we observe that the 
region v>A is never reached by the system, as showed in the inset of Fig. 1; this region is characterized by a'(v)<0. 
We see that: i) for fc=l, where there is no phase transition, the function a(v) is analytic; ii) for k=2, when we 
observe a second order PT, the first derivative of o~(v) is discontinuous at v c =v(e c )=A, and its second derivative is 
negative around the singular point. Hi) for fc>3 the first derivative of cr(v) is also discontinuous at the transition point 
v c =A, but its second derivative is positive around v c . In this case a first order transition takes place. Therefore the 
investigation of the potential energy topology, via <j(v), allows us to establish the location and the order of the PTs, 
without introducing any statistical measure. 

The previous results allows us to conjecture that there is a relation between the thermodynamic entropy of the 
system and <j(v). In fact, the presence of a first order transition with a discontinuity in the energy is generally related 
to a region of negative specific heat, i. e. of positive second derivative of the entropy. Thus, it seems that at least 
around the transition point the thermodynamic entropy and o~{v(e)) are closely related, in the sense that the jump in 
the second derivative of s(e) is determined by the jump in the second derivative of o~(v(e)). Then it should be possible 
to write 




(25) 




(26) 



s(e) ~ <r(v(e)) + 71(e) 



(27) 
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where IZ(e) is analytic (or, at least, C 2 ) around the transition point. We hope that future work will address this point 
in a more quantitative way. 

In conclusion, we have shown that the presence of a phase transition and its order are signaled in a very clear way 
in the topology changes of the potential energy hyper-surfaces. In our model, the Euler characteristic of the region 
of phase space in which V(q)<Nv shows a discontinuity in its first derivative with respect to v if and only if a, phase 
transition takes place at the same energy value. From the concavity of the Euler characteristic as a function of v 
around the singular point we can deduce the order of the transition, that is first order if d 2 x(v)/dv 2 >0 and second 
order if d 2 \/ 'dv 2 <0. From the last observation we conjecture a deep relation between the thermodynamic entropy 
and some topological property of the potential energy landscape, at least around the transition point. 

We acknowledge support from INFM and MURST COFIN2000. 
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